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The field theoretic renormalization group and the operator product expansion are applied to the 
model of passive vector (magnetic) field advected by a random turbulent velocity field. The latter is 
governed by the Navier-Stokes equation for compressible fluid, subject to external random force with 
the covariance oc S(t — t')k'^~‘^~^, where d is the dimension of space and y is an arbitrary exponent. 
From physics viewpoints, the model describes magnetohydrodynamic turbulence in the so-called 
kinematic approximation, where the effects of the magnetic field on the dynamics of the fluid are 
neglected. The original stochastic problem is reformulated as a multiplicatively renormalizable field 
theoretic model; the corresponding renormalization group equations possess an infrared attractive 
fixed point. It is shown that various correlation functions of the magnetic field and its powers 
demonstrate anomalous scaling behavior in the inertial-convective range already for small values 
of y. The corresponding anomalous exponents, identified with scaling (critical) dimensions of certain 
composite fields (“operators” in the quantum-field terminology), can be systematically calculated as 
series in y. The practical calculation is performed in the leading one-loop approximation, including 
exponents in anisotropic contributions. It should be emphasized that, in contrast to Gaussian 
ensembles with finite correlation time, the model and the perturbation theory presented here are 
manifestly Galilean covariant. 
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I. INTRODUCTION 

Much attention has been attracted to the problem of 
intermittency and anomalous scaling in developed mag¬ 
netohydrodynamic (MHD) turbulence; see, e.g., [lj|-[^ 
and references therein. It has long been known that 
in the so-called Alfvenic regime, the MHD turbulence 
demonstrates the behavior, similar to that of the usual 
fully developed fluid turbulence: cascade of energy from 
the infrared (IR) range towards smaller scales, where the 
dissipation effects dominate, and self-similar (scaling) be¬ 
havior of the energy spectra in the intermediate (inertial- 
convective) range. Moreover, strongly non-Gaussian (in¬ 
termittent) character of the fluctuations in the MHD tur¬ 
bulence is much strongly pronounced than in ordinary 
turbulent fluids or in the passive scalar problem. 

The solar wind, a conducting fluid expanding into the 
interplanetary space, covers a wide range of spatial and 
temporal scales and thus provides a kind of “laboratory” 
in which various models of the MHD turbulence can 
be tested 0“[il|- In solar flares, highly energetic and 
anisotropic large-scale events (with the magnetic fields as 
intense as 500 Gauss) coexist with small-scale stochastic 
fluctuations and coherent structures, finally responsible 
for the dissipation. Thus, modelling the way how the en¬ 
ergy is distributed, conveyed along the scales and even- 
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tually dissipated is a difficult task. 

The intermittency strongly modifies the IR behavior of 
the higher-order correlation functions, leading to anoma¬ 
lous scalin g w ith infinite sets of independent anomalous 
exponents Q. 

A simplified description of the situation was proposed 
in the large-scale field = UiB^ dominates the 
dynamics in the distinguished direction specified by a 
unit constant vector n = {rii}, while the fluctuations 
in the perpendicular plane are described as nearly two- 
dimensional. This picture allows for precise numerical 
simulations, which show that turbulent fluctuations or¬ 
ganize in rare coherent structures separated by narrow 
current sheets. On the other hand, the satellite obser¬ 
vations and numerical simulations sa suggest that 
the scaling behavior in the solar wind is closer to the 
anomalous scaling in the three-dimensional fully devel¬ 
oped hydrodynamic turbulence, rather than to simple 
Iroshnikov-Kraichnan scaling a a suggested by the two- 
dimensional picture with the inverse energy cascade. 

Thus, further analysis of more realistic three- 
dimensional models is welcome. 

In a number of papers, the problem was studied within 
the framework of the kinematic approximation, in which 
the magnetic field is passive in the sense that it does not 
affect the dynamics of the velocity field [i 3-[23. This 
approximation seems reasonable if the gradients of the 
magnetic fields are not too large. What is more, the 
renormalization group analysis of [s^ suggests that such 
a “kinematic regime” can indeed describe the possible IR 
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behavior of the full-scale problem. It is then possible to 
model the velocity field “by hands,” that is, by simple 
statistical ensembles with prescribed properties. Most 
popular is the Kazantsev-Kraichnan ensemble [3l|, [s^: 
the random velocity field is Gaussian, white in time, and 
has a power-law spectrum. 

Numerous analytical and numerical results were de¬ 
rived for the scalar and vector fields, advected by the 
Kazantsev-Kraichnan “flow,” see [s^ for the review and 
references. The main results concerning anomalous scal¬ 
ing for the magnetic case can be summarized as follows 

[il-lll: 

(i) Anomalous scaling is present and appears already 
for the pair correlation function. 

(ii) In the presence of large-scale anisotropy (brought 
about, e.g., by the constant background field 5°), the 
anomalous exponents for a given correlation function 
demonstrate a kind of hierarchy: in the expansion of cor¬ 
relation functions in the spherical harmonics Yim, the 
corresponding exponents increase with Z, the degree of 
anisotropy. Thus, for the even-order functions, the lead¬ 
ing terms of the inertial-range behavior are given by the 
isotropic contribution {I = 0). This gives quantitative 
support for Kolmogorov’s hypothesis of the local isotropy 
restoration. 

(iii) Nevertheless, the anisotropy survives at small 
scales and manifests itself in odd-order correlation func¬ 
tions, or in dimensionless ratios involving such functions 
(like the skewness factor). 

An important advantage of the Kazantsev-Kraichnan 
ensemble is the possibility to easily model anisotropy 
and compressibility. Importance of compressibility for 
the MHD turbulence was realized already in the classical 
study of [^. Within the framework of the Kazantsev- 
Kraichnan ensemble, effects of compressibility were stud¬ 
ied in [3, It was shown that: 

(iv) The anomalous exponents depend on the degree 
of compressibility. When it grows, the hierarchy of 
anisotropic contributions becomes less pronounced and 
the persistence of anisotropy in the depth of the inertial 
interval becomes more remarkable. 

Of course, generalization of this analysis to more real¬ 
istic velocity dynamics is necessary: some of the afore¬ 
mentioned results can be artefacts of the oversimplified 
statistics. 

It is possible to directly generalize the Kazantsev- 
Kraichnan ensemble to the case of finite correlation time; 
see, e.g., for the passive scalar and for the 

passive vector fields. However, such “synthetic” models 
with non-vanishing correlation time suffer from the lack 
of Galilean symmetry, which may lead to “interesting 
pathologies,” quoting the authors of [13 ■ One of such a 
pathology manifests itself as ultraviolet (UV) divergence 
in the vertex [35l| - which in more realistic models is for¬ 
bidden by Galilean invariance. 

Thus, it is desirable to describe the advecting veloc¬ 
ity field by the Navier-Stokes equations with a random 
stirring force and to work within Galilean covariant for¬ 


malism. For the incompressible case, the analysis of the 
passive vector field was accomplished in [29j |- 

In this paper, we study the anomalous scaling in the 
kinematic MHD problem and model the velocity dynam¬ 
ics by the non-Gaussian velocity field with finite correla¬ 
tion time, governed by the stochastic Navier-Stokes (NS) 
equation. We apply to the problem the approach based 
on the field theoretic renormalization group (RG) and 
the operator product expansion (OPE), earlier applied 
to the passive scalar problem . In that approach, 

the anomalous exponents are identified with the critical 
dimensions of certain Galilean-invariant composite fields 
(“operators”). It can be directly generalized to the cases 
of finite correlation time, presence of anisotropy, non- 
Gaussianity and so on. Passive advection of vector fields 
(and hence kinematic MHD problems) with various veloc¬ 
ity ensembles were studied earlier within the RG-I-OPE 
approach in [^-[^. 

A general overview of the RG-I-OPE approach to the 
problem of anomalous scaling and more references can 
be found in [ 43 . Detailed exposition of earlier work on 
the field theoretic RG approach to stochastic models of 
turbulence on the whole is presented in [4l|- 

However, analysis of the compressible fluid on the base 
of stochastic NS e quat ion appears a difficult task in it¬ 
self; see, e.g., [^-[^. In spite of some discrepancies, 
all of those studies support the existence of a “strongly 
compressible” scaling regime, different from the usual in¬ 
compressible one. 

In the present paper, we adopt the approach of [T7I| . 
where, with the price of some natural approximations, 
the stochastic NS equation for a compressible fluid was 
reformulated as a multiplicatively renormalizable field 
theoretic model. Then the standard field theoretic RG 
was applied to the problem, and the resulting stationary 
scaling regime was associated with the IR attractive fixed 
point of the corresponding RG equations. 

Recently, that ensemble was employed to study, within 
the RG-I-OPE framework, the problem of passive scalar 
advection in a turbulent compressible fluid [^ . In spite 
of close resemblance with the case of Kraichnan’s model, 
some of the results appeared somewhat different. The 
present paper continues the study of in connection 
with the MHD turbulence. Eor this reason, we will only 
briefly discuss the points, common to the scalar and vec¬ 
tor problems, refer the reader to the papers [13,113 when¬ 
ever possible, and focus on the points specific of the vec¬ 
tor case. 

The plan of the paper is the following. In section [IT] 
we give the detailed description of the model: the veloc¬ 
ity ensemble, the stochastic MHD equation and the field 
theoretic formulation. In section Ell we discuss canon¬ 
ical dimensions and renormalization of the field theo¬ 
retic model, demonstrate its multiplicative renormaliz- 
ability and calculate (in the leading one-loop approxi¬ 
mation) the corresponding renormalization constant. In 
section lTVl we derive the corresponding RG equations and 
show that they possess the only IR attractive fixed point 
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in the physical region of parameters. This fact implies 
the scaling behavior in the IR range (long times, large 
distances); the corresponding critical dimensions of the 
basic fields and parameters are presented. In section El 
we calculate the critical dimensions of the tensor com¬ 
posite fields (operators), constructed solely of the basic 
scalar fields; these will play the crucial role in the fol¬ 
lowing. In section |VT] we employ the OPE to derive the 
inertial-range asymptotic behavior of various correlation 
functions. Section IVlIl is reserved for discussion and the 
conclusion. 

II. THE MODEL 
A. The Velocity Ensemble 

Following [U , we describe the stochastic dynamics 
of a compressible fluid by the set of two equations: 

VtUi = - didk\vk + iJ.odidkVk -di4>+fi (2.1) 

Vt^ = -cld,v^, (2.2) 

which are derived from the momentum balance equation 
and the continuity equation with two assumptions: 
the kinematic viscosity coefficients vq and are assumed 
to be constant, that is, independent oi x = {t,x}, and 
the equation of state is taken in the simplest form of the 
linear relation {p — p) = Cq{p — p) between the devia¬ 
tions of the pressure p(x) and the density p(x) from their 
mean values; then the constant Cq has the meaning of the 
(adiabatic) speed of sound. 

In Eqs. (12.11) . (12.2L v = {vi{x)} is the velocity field 
and, instead of the density, we use the scalar field defined 
as (j){x) = Cq ln(p(a;)/p). Furthermore, 

= dt+ Vkdk (2.3) 

is the Lagrangean (Galilean covariant) derivative, dt = 
d/dt, di = djdxi, and 9^ = dA is the Laplace operator. 
The problem is studied in the d-dimensional (for gener¬ 
ality) space X = {xi}, i = 1 .. .d, and the summations 
over the repeated Latin indices are always implied. 

In the Navier-Stokes equation fi is the density 

of the external force (per unit mass), which mimics the 
energy input into the system from the large-scale stir¬ 
ring. In order to apply the standard perturbative RG to 
the problem, it is taken to be Gaussian with zero mean, 
not correlated in time (this is dictated by the Galilean 
symmetry), with the given covariance 

{h{x)fj{x')) = 6{t - t') f D{.{k) expjik ■ x}, 

Jk>m 

(2.4) 

with 

D{^(k) = Do {/^|(k) + aPii(k)} . (2.5) 

Here (k) = 6ij — kikj /fc^ and (k) = kikj / k"^ are the 
transverse and the longitudinal projectors, respectively. 


k = |k| is the wave number (momentum), Dq and a are 
positive amplitudes. The parameter go = Dojv^ plays 
the part of the coupling constant (expansion parameter 
in the perturbation theory); the relation go de¬ 

fines the typical UV momentum scale. The parameter 
m ~ reciprocal of the integral turbulence scale, 

provides IR regularization; its precise form is unessen¬ 
tial and the sharp cut-off is merely the simplest choice 
for calculational reasons. The exponent 0 < y < 4 plays 
the part analogous to e = 4 — d in the RG theory of 
critical state : it provides UV regularization (so 

that the UV divergences have the form of the poles in 
y) and various scaling dimensions are calculated as series 
in y. The most realistic (physical) value is given by the 
limit ?/ —>■ 4: then the function (1^ can be viewed as 
a power-like representation of the function i5(k) and it 
corresponds to the idealized picture of the energy input 
from infinitely large scales. 

As already mentioned, more detailed justification and 
discussion of the compressible model (I2.1I) - (I2.5I) is given 
in PS. 

B. The MHD Equation 

In the presence of a constant background field = 
B^rii with a certain constant unit vector n = {rii}, the 
dynamic equation for the fluctuating part 6i = Oi{t, x) of 
the full magnetic field Bi = B^{ni -\- 61 ) has the form 

dtOi + dk{vk9r - OkVi) = Kod^Oi + UkdkV^, ( 2 . 6 ) 

where kq = c^I^ko is the magnetic diffusion coefficient. 
Equation (12.61) follows from the Maxwell equations ne¬ 
glecting the displacement current and the simplest form 
of Ohm’s law for a moving medium j = tr (E -|- [u, B]), 

where a is the conductivity and c/ is the speed of light; 
see, e.g., [soj - 

The last term on the right-hand side of (1^ main¬ 
tains the steady state of the system and acts as the 
source of the anisotropy; in principle, it can be replaced 
with an artificial Gaussian noise with appropriate statis¬ 
tics. In the real problem, the field v satisfies the Navier- 
Stokes equation with the additional Lorentz force term 
^ (B X curlB). In our kinematic approximation the back 
reaction of the magnetic field on the velocity dynamics 
is neglected and the latter is described by the stochastic 
problem (I2.1I) - (I2.5I) without the Lorentz term. 

C. The field theoretic formulation 

It is well known that any stochastic problem of the 
type (l^ - ((23)) can be reformulated, in a standard fash¬ 
ion, as a certain field theoretic model; see, e.g., [5l|,P- 
This means that various correlation and response func¬ 
tions of the original stochastic problem can be repre¬ 
sented as functional integrals over the doubled set of 
fields 4) with the weight exp5(4)), where iS($) is the 
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so-called De Dominicis-Janssen action functional. The 
action functional 5„(<i>) for the problem (I2.1I) - (I2.5I) with 
$ = {v', (j)', V, (j)} looks too cumbersome, and we do not 
reproduce it here, as well as the elements of the corre¬ 
sponding Feynman diagrammatic techniques (bare prop¬ 
agators and vertices); they can be found in [^ |^. Be¬ 
low we only will need the velocity-velocity propagator at 
Co = 0; in the frequency-momentum (w-k) representa¬ 
tion it has the form: 


{viVj)o = Do 


up- -I- 


1 

uP + uliu'^k^ J ■ 


(2.7) 


The full-scale stochastic problem (lOl-dTHI) corre¬ 
sponds to the action functional 

5($) =5,(-u',0',-u,0)+5e(0',e,-u), (2.8) 


where 


^0 = 9'^{-dt6i-dkivkOi-dkVi) + Kod'^9i + nkdkVi} (2.9) 

is the De Dominicis-Janssen action for the problem (12.61) 
at fixed v. It brings about the new propagator 


mu = 


( 2 . 10 ) 


—iw + wovok^ 
and the new vertex Viji9[9jVi with the vertex factor 

Viji{k) = \{5ijki - Stikj). ( 2 . 11 ) 


dp = dp + 2d^, dehned such that all the viscosity or 
diffusivity coefficients are dimensionless; see [4l|, [52|. All 
the canonical dimensions in our model (EDl-dSll) are 
identical to their counterparts in the scalar case, and we 
refer the reader to the Table 1 in [j^ . 

The formal index of UV divergence of a certain 1- 
irreducible Green function T is given by its total canoni¬ 
cal dimension : 

6r = d 2 — 'y ^ N^d^, (3.1) 

3> 

where are the numbers of the fields entering into the 
function T, are their total canonical dimensions. Su¬ 
perficial UV divergences, whose removal requires coun¬ 
terterms, can be present only in the functions T with a 
non-negative integer Jr- The counterterm is a polyno¬ 
mial in frequencies and momenta of degree Jr, with the 
convention that uj ^ k'^. 

The dimensional analysis (“power counting”) should 
be augmented by certain additional considerations: 

(i) All the 1-irreducible Green functions without the 
auxiliary (“primed”) fields vanish identically and thus 
require no counterterms. 

(ii) If a number of external momenta occurs as an over¬ 
all factor in all the diagrams of a certain Green function, 
the real index of divergence Jp is smaller than Jr by the 
corresponding number of unities. In the model Sy the 
field (f> enters the corresponding vertex only in the form 
of a spatial derivative, which reduces the real index of 
divergence: 


A few remarks are in order here. First, the derivative 
at the vertex in (12.911 can be moved onto the auxiliary 
field 6' using the integration by parts; thus k in (12. lip is 
the momentum argument of O'. Second, the vertex factor 
satisfies the transversality condition 

fc,Uy,(k)=0, (2.12) 

that follows from its explicit form (12.lip . It remains to 
note that another new mixed propagator {9v)o will not 
appear in relevant diagrams. 

III. UV DIVERGENCES AND THE 
RENORMALIZATION 

A. Canonical dimensions, connterterms, and 
renormalizability 


J^ = Jp - N^. (3.2) 

(iii) The Galilean invariance of the model requires that 
the counterterms be also invariant. In particular, this 
means that the covariant derivative (|2.3I) appears in the 
counterterms as a whole. 

These considerations forbid superficial UV divergences 
in certain Green functions, allowed by dimensional anal¬ 
ysis, and hence reduce the number of the counterterms 
needed for the renormalization of the model. 

The analysis of the field theoretic model with the ac¬ 
tion Sy in (12.8p . performed in (see also [1^), has 
shown that it is multiplicatively renormalizable (after a 
simple natural extension). This means that all the UV 
divergences can be removed from the Green functions by 
the renormalization of the fields (f> m Z^cj), (j)' m 
and of the parameters: 


The analysis of UV divergences is based on the anal¬ 
ysis of canonical dimensions; see, e.g., [sllH^. Dynam¬ 
ical models of the type (IO)-(l^ have two scales: the 
time scale T and the length scale L. The canonical di¬ 
mension of any quantity F (a field or a parameter) is 
described by two numbers, the momentum dimension 
dp and the frequency dimension dp, defined such that 
[F] ^ . In the renormalization theory, the 

central role is played by the total canonical dimension 


go — 90^Zg, voZuj cq — cZ^ , (3.3) 

and so on. Here the renormalization constants Zi absorb 
all the UV divergences, so that the Green functions are 
UV finite (that is, finite at y = 0) when expressed in 
terms of the renormalized parameters g, u, and so on; the 
reference scale (or the “renormalization mass”) g, is an 
additional free parameter of the renormalized theory. No 
renormalization of the fields v',v and of the parameters 
TO, a is needed. 
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The inclusion of the new contribution Sg in the full 
model brings about the only new UV divergence in the 1- 
irreducible function { 0 ' 0 )i-ir with the counterterm O'd^O. 
Two points are important here: 

(iv) From the linerarity of the original stochatic model 
in the field 9 it follows that Ngi —Ng is a. non-negative in¬ 
teger for any nontrivial 1-irreducible Green function: no 
other Feynman diagram can be drawn. This fact forbids 
the superficial divergences in all the 1-irreducible func¬ 
tions {9'0 ... 6 )i-ir, except for the first one, and thus pre¬ 
vents our model from being non-renormalizable, despite 
the fact that the magnetic field has a negative canonical 
dimension. 

(v) For the full model (|2.8I) . the items (ii) and (iii) 
require some additional discussion. The derivative at the 
vertex in Sg can be moved, using the integration by parts, 
onto the field O'. Thus, the real index of divergence is 
reduced according to the item (ii) above, and 9' enters the 
countertems only as a spatial derivative. The expression 
(13.21) has to be replaced with 


5'^ = 5^-N^-Ng,. (3.4) 


Thus, the counterterm O'dtO is forbidden, and so is 
9'{vidi)9 due to the Galilean symmetry (iii). 

The only remaining counterterm 9'd^9 is naturally re¬ 
produced by the multiplicative renormalization of the 
magnetic diffusion coefficient: kq = kZ,^. No renormal¬ 
ization of the fields 9', 9 is needed. 

The renormalized analog of the action functional (12.81) 
has the form 


5^(<l>)=5«(4>)+5«($), (3.5) 


where 5'^($) is the renormalized analog of the action 
5($), given in [13,13, and 


^g^ = (^i{-dt9i - dk{vk9i - 9kVi) + nZ^d'^Oi + nkdkVi} 

(3.6) 

is the renormalized part of the full action that describes 
interaction with the magnetic field. 

I 


B. Leading-order calculation of the 
renormalization constant 


We performed the practical calculation of the renor¬ 
malization constant in the leading one-loop approxi¬ 
mation, which is consistent with the accuracy of the cal¬ 
culation for the NS problem dUD, dH made in (^ . 
Although this calculation is rather simple and similar 
to that for the Gaussian velocity ensemble 2^, we will 
briefly discuss it for the sake of completeness and in order 
to stress some peculiarities. 

The constant is found from the requirement that the 
1-irreducible Green function {9'9)i-iY be UV finite (that 
is, finite at y ^ 0) when expressed in renormalized pa¬ 
rameters. In the frequency-momentum representation it 
has the form: 


(6l(6l2)l-ir(U,p) = {-Kop^-f iU} Pf^(p)-b Ei2(U,p), 

(3.7) 

where S 12 is the “self-energy operator” given by infinite 
sum of 1-irreducible Feynman diagrams and p = |p|. Be¬ 
cause of the large number of tensor indices involved in 
our expressions, we use numbers (instead of latin letters) 
to denote them, with the standard convention on the 
summation over repeated indices. 

The only one-loop self-energy diagram looks as follows: 



Here the wavy line denotes the bare propagator {vv)o, the 
solid line with a slash denotes the bare propagator {99')q 
from (I2.10|) . the slashed end corresponding to the field 
9'. The dots with three attached fields 9', 9, v denote 
the vertex (12.111) . 

In this approximation, the renormalization constant 
in the bare term of (13.71) is taken to the first order in g, 
while in the diagram (j3.8l) all Z’s are simply replaced with 
unities. Furthermore, we only need to know the divergent 
part of (1^ . which is quadratic in p (see the preceding 
subsection). Thus, we can put H = 0 in (j3.7l) and retain 
only quadratic terms in the expansion of Ei 2 (H = 0,p) 
in p. Like for the original NS model, its divergent part 
is independent of cq ^ c and can be calculated directly 
at c = 0; see the discussion in [3l. Thus, we can use the 
expression (E31) for {vv)o- 

Then the analytic expression for (|3.8I) takes on the 
form: 


Ei 2 (H — 0 , p) — Dq 


duj 

27r 


' k>m 


dk 

(27ry 


Vl43(p)G625(P + k) 


uj'^ + 


g-PiUk) \ p4i(k + p) , , 

Lu^ + j —iw-b -b pp ’ 


where k = |k|. The simplifying replacement 64 ^ dition (12.121) and the contraction with V 625 . 

can immediately be made owing to the transversality con- Integrations over the frequency are easily performed. 
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for example: 


duj 1 1 1 

27r —iw + + kp u^v'^k'^ 2ui''^kP'{uk'^ -\- w\tp + \< l \'^) 


(3.10) 


The numerators in the integrand of (13.91) contain the 
terms quadratic and linear in p. For the first ones, one 
can immediately set p = 0 in (13.101) . while for the second 
ones, one should expand (13.101) up to the linear term in 
p, for example. 


the resulting integrals are reduced to the scalar integral 

= (3.12) 

where 


_1_^ 1 L_ 2w (pk) \ 

ufc^ + w|p + kp {u + w)k^\ (u + w) j 

With the aid of the formulas 

Jdkhfik) = 0, Jdk^fik) = ^-^ J dkfik), 

J dkfik), 

(3.11) 

where f{k) is any function depending only on fc = |k|, all 

I 


Sd = 27r'^/Vr(d/2) (3.13) 

is the surface area of the unit sphere in the d-dimensional 
space and r( - • •) is Euler’s Gamma function. 

The final result contains two types of terms, propor¬ 
tional to 1^12 (p) and p| 2 (p)i respectively. Owing to the 
transversality of the fields 0, O', the latter ones should be 
discarded. (This would happen automatically if we in¬ 
cluded the corresponding projector into the vertex (12.111) , 
but we did not do that for brevity). Practically, it is more 
convenient to collect only terms proportional to 612 
and drop all the terms proportional to P 1 P 2 in the course 
of calculation. It remains to express the amplitude Dq in 
(13.9p in terms of renormalized variables: Dq = . 

Then the final result reads: 


dk 


k^kgkikp r/1\ _ digSip -f SiiSgp S'lpSgi 


k^ 


-fik) = 


d{d + 2 ) 


Ei2(n = 0,p) 


-iyp^P^2iP) 


9 

2 dy 



jd-l) 

1 -I- u> 


a{u — w) 
u{u + w)^ 


(3.14) 


Here we passed to the new coupling constant 

g = gSd/i2TT)'‘, (3.15) 

with Sd from (13.131) . 

Then in the MS scheme the renormalization constant 
Zk that cancels the pole of the expression (13.141) in the 
renormalized analog of the function (13.71) (that is, with 
the replacement kq kZ^ in the bare term) has the 
form: 


by the same velocity ensemble; see Eq. (3.24) in (4^ . 
Similar coincidence between the passive scalar and mag¬ 
netic fields in the kinematic approximation was earlier 
observed for the incompressible case (see, e.g., 1^) : 
sometimes it extends to the two-loop approximation [28l| . 

IV. RG EQUATIONS, FIXED POINT AND THE 
CRITICAL DIMENSIONS 


g r jd-l) a(u - w) ) 

2 dwy \ (1 -I- w) u{u -I- w)^ J ’ 


(3.16) 


while the corresponding anomalous dimension is 


g J (d — 1) a{u — ui) 1 

2 dw \ (1 -I- w) u{u + wy J ’ 


(3.17) 


with the corrections of the order g^ and higher. 

It is interesting to note that the expression (13.171) liter¬ 
ally coincides with its analog for the scalar fields advected 


Here we only briefly discuss the derivation of the IR 
scaling behavior from the RG equations in our model; it 
is nearly identical to the scalar case, discussed in [4^ in 
great detail. 

Multiplicative renormalizability of the field theoretic 
model (12.81) allows one to derive, in a standard way, dif¬ 
ferential RG equations for the renormalized Green func¬ 
tions 


G{e,p,...) = ($...$)fi. 

Here e = {g,iy,u,v,w,c,m,a} is the full set of renor- 
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malized parameters, /r is the reference momentum scale 
and the ellipsis stands for the other arguments (times 
or frequencies and coordinates or momenta). For con¬ 
venience, we introduced here three dimensionless ratios: 
uo = and vg = Xo/vg are related to the viscos¬ 

ity and diffusivity coefficients of the (properly extended) 
model ( 1 ^ . ( 1 ^ . while wg = Kg/ug is related to the 
magnetic diffusivity coefficient; u,v,w are their renor¬ 
malized analogs. 

The RG equation expresses the invariance of the renor¬ 
malized Green function with respect to changing of the 
reference scale /r, when the bare parameters eg are kept 
fixed: 


+ G(e,M,...)=0. (4.1) 

Here and below we denote Vx = xdx for any variable x 
and is the operation = ^9^ at hxed eg. In terms 
of the renormalized variables, it takes the form 

T’/i = + Pg9g + j3udu + Pvdv + Pwdw — Iv'Dv — XcDc- 

(4.2) 

The anomalous dimension xf of a certain quantity F (a 
field or a parameter) is defined by the relation 


7 F='D^lnZF, (4.3) 

and the /3 functions for the dimensionless parameters 
(“coupling constants”) are 

='^ti9 = 9[-y -Ig]^ 

Pu = Vfj.u = -U7„, (4.4) 

and similarly for I3y, j5w Here the second equalities result 
from the definitions and the relations of the type (j3.3l) . 

Note that from the definition of wg it follows that Zi^ = 
Z^Z.un so that 


fSw = w[xv - lw\- (4.5) 

The possible types of IR asymptotic behavior are as¬ 
sociated with IR attractive fixed points of the RG equa¬ 
tions. The coordinates g* = {gi*} of the fixed points are 
found from the equations 

ft(5,)=0, (4.6) 

where g = {gi} is the full set of coupling constants and 
Pi = F^gi are their P functions. The character of a fixed 
point is determined by the matrix 

Vlij = dPi/dgj\g=g^. (4.7) 

For the IR fixed points the matrix H is positive (that is, 
positive are the real parts of all its eigenvalues). 

The analysis performed in (see also [d^) on the 
base of the leading-order (one-loop) approximation has 


shown that the RG equations of the model corre¬ 
sponding to the stochastic NS problem (EH), EH, pos¬ 
sess the only IR attractive fixed point in the physical 
region of parameters {g,u,v > 0 ): 

5 * =-^^^^^ + 0 ( 2 /^), u* = l-bO(?/), v., = l + 0{y). 

(4.8) 

From a certain exact relation between the renormal¬ 
ization constants [d^, the exact result 

7 * = y/3 (4.9) 

follows (no corrections of the order y^ and higher). Here 
and below, 7 * denotes the value of the anomalous dimen¬ 
sion 7 i at the fixed point. 

Now we substitute the one-loop expressions (I3.17L 
(EH) and the exact result (EH) into Eq. (14.51) . Then 
the equation /3 k, = 0 yields, after simple algebra, the 
equation 

{w — l)[(c? — l)(w -I- l)(w -I- 2) -I- 2 q;] = 0, (4-10) 

which has the only positive solution w* = 1 , with possible 
corrections of order 0 {y) and higher. 

Since the functions Pg^u,v do not depend on w, the new 
eigenvalue of the matrix (14.71) coincides with the diagonal 
element 


dp.u,ldw\g=g, = y {3((i - 1) -f a} /&{d - 1) > 0; 


thus the fixed point ()4.8I) and = 1 of the full model 
remains IR attractive. 

Existence of an IR attractive hxed point in the physi¬ 
cal region of the parameters implies existence of scaling 
behavior in the IR range. The critical dimension of some 
quantity F (a held or a parameter) is given by the rela¬ 
tion (see [ 41 . [52|1 


= 4 + AK,d‘^ + 7F, Ak, = 2 - 7 : = 2-y/3. (4.11) 


Here dPp and are the canonical dimensions of E, Xf 
is the value of the anomalous dimension xf at the hxed 
point, and A^, is the critical dimension of the frequency. 

The critical dimensions of the helds and parameters 
of the model described by the action Sy from Eq. (12.81) 
are presented in [i^; see also [i^. In addition, our full 
model involves two more critical dimensions: 


Ae = -l-|-y/ 6 , Ag/= d-I-1 - 2 // 6 . (4.12) 

These expressions are exact because the helds 9 and 9' 
are not renormalized. 


V. COMPOSITE FIELDS AND THEIR 
DIMENSIONS 

An important role in the following will be played by 
certain composite helds (“composite operators” in the 






quantum-field terminology). In general, a local compos¬ 
ite operator is a monomial or polynomial built of the 
primary fields $(a;) and their finite-order derivatives at 
a single space-time point x. In the Green functions with 
such objects, new UV divergences arise due to coinci¬ 
dence of the field arguments. They should be eliminated 
by additional renormalization procedure. As a rule, op¬ 
erators mix in renormalization: renormalized operators 
(whose Green functions are UV finite) are given by fi¬ 
nite sums of the original monomials. However, in the 
following only a simpler situation will be encountered, 
when the original operator F{x) and its renormalized 
analog F^{x) are related by multiplicative renormaliza¬ 
tion F{x) = ZfF^{x) with a single renormalization con¬ 
stant Zp- Then the critical dimension Af of the oper¬ 
ator F is given by the same expression (14. Ill) and, in 
general, differs from the naive sum of the dimensions of 
the fields and derivatives that compose the operator. 

We will focus on the irreducible tensor operators built 
solely of the fields 9. They have the forms 

Fni = Oi^ (x)--- 9i, {x) {0i{x)9i{x)y + (5.1) 

where I < n is the number of free tensor indices and 
n = I + 2s is the total number of the fields 9 entering into 
the operator; the tensor indices and the argument x of 
the symbol Fni are omitted. The ellipsis stands for the 
appropriate subtractions involving the Kronecker delta 
symbols, which ensure that the resulting expressions are 
traceless with respect to contraction of any given pair of 
indices, for example, 9i9j — Sij{9k9k/d) and so on. 

The total canonical dimension of any 1-irreducible 
Green function T with one operator F(x) and arbitrary 
number of primary fields (the formal index of UV diver¬ 


gence) is given by 

(5r = dp ~ 'y ^ N^d^, (5.2) 

$ 

where are the numbers of the fields entering into T, 
d^ are their total canonical dimensions, and dp is the 
canonical dimension of the operator. Superficial UV di¬ 
vergences can be present only in the functions T with 
a non-negative integer 6 r- For the operators (15.ip from 
Table I in we find dp = —n. The linearity of the 
equation (12.61) in the field 9 imposes the restriction that 
Ne in (j5.2p cannot exceed the number of the helds 9 in 
the operator F. The direct analysis shows that super¬ 
ficial UV divergences (dr > 0) for Fni can be present 
only in the 1-irreducible functions with Ngi = Ny = 0 
and Ng = n; they are all logarithmic: dr = 0. The 
simple inspection shows that the mixed propagator {9v)o 
does not appear in the relevant Feynman diagrams; in 
other words, the last term in the right-hand side of the 
equation (12.61) is unimportant here. Without that term, 
the model becomes 0{d) invariant. In turn, this means 
that irreducible operators with different values of I cannot 
mix with each other. We finally conclude that the oper¬ 
ators (EH) renormalize multiplicatively: Fni = ZniF^i 
and turn to the one-loop calculation of the renormaliza¬ 
tion constant Zni and of the critical dimension of the 
operator EH, which will be denoted as A„/. 

Let r(a:;0) be the generating functional of the 1- 
irreducible Green functions with one composite operator 
F{x) = Fni and any number of fields 9. Here x = {t,x} 
is the argument of the operator and 9 is the functional 
argument, the “classical analog” of the random field 9. 
We are interested in the n-th term of the expansion of 
F(a:; 9) in 9, which we denote F„(a:; 9). It can be written 
as follows: 


r„(x;0) 



dXn 9{xi) ■ ■ ■ 9{Xn) {F{x)9{xi) ■ ■ ■ 9{Xn))l-{v 


(5.3) 


In the one-loop approximation the function (15.3p is rep¬ 
resented diagramatically as follows: 



The first term is the tree (loopless) approximation, and 
the thick dot with the two attached lines in the diagram 
denotes the operator vertex, to be specified later. 

The renormalization constant Zni for the operator 
EH is found from the requirement that the renormalized 
analog F^ = Zi^^Fn of the function (15.3|) be UV finite in 
terms of renormalized parameters. 

For practical calculations, it is convenient to con¬ 


tract the tensors EH with an arbitrary constant vector 
A = {Ai}. The resulting scalar operator has the form 

=(A,0,)'(0,0,)* + ..., (5.5) 

where the subtractions, denoted by the ellipsis, necessar¬ 
ily involve the factors of A^ = AAf 

Within our accuracy, it is sufficient to replace all the 
renormalization constants in the diagram with unities, so 
that Dq gv^, uq ^ u and so on. Furthermore, we 
are eventually interested in the fixed-point value of the 
anomalous dimension, so that we can set rt = re = 1 in 
the following. Since the diagram is logarithmically diver¬ 
gent, we can set all the external frequences and momenta 
equal to zero. 

Like for the calculation of the self-energy diagram in 
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sec. IIII B[ here and below we use numbers (instead of 
Latin letters) to denote the tensor indices. Then the 
diagram in (15.411 can be represented as follows: 


for the transverse contribution in (|2.7I) and 
f dk 1 


lk>m ( 2 ^)'" 


PfV(k)P 2 i(k) 


(5.11) 


Vi2{0) (71278 ^7^8, 


(5.6) 


where V12 {&) is the operator vertex (denoted by the thick 
dot in the diagram and specified below), the fields 979 s 
(denoted by wavy tails) are attached to the lower vertices 
(12.111) (small dots) and (7i278 is the “core” of the diagram. 
It has the form 


for the longitudinal one. With the aid of the relations 
(13. Ill) , all these integrals are reduced to the scalar inte¬ 
gral (13.121) . 

Combining all these contributions and contracting the 
result with the fields 0708 gives for (15.6p the following 
expression: 


-Dn 


■ Vi2 {Ti 2 -I- CXL12} . 


(7i278 — 


dui 


dk gv'^^yRi27s 


27 r Jk>m ( 27 r)‘^ (w^ -h 


(5.7) 


with the tensorial factor stemming from the vertices 
(12.111) and the projectors of the propagators: 


-Ri 278 = Pii(k)P2i(k) |P5'6(k) + aPiUk)} 1475(k)F486(k). 

(5.8) 

Intergation over the frequency is easily performed: 


where 


and 


Ad{d + 2)v^ 

Ti2 = {d l)di20^ ~ 20102 


(5.12) 


P12 — < 5 i 20 ^ + { d ^ — 2)0102. 
Now let us turn to the vertex factor 

V 12 = 


dui 


2 tt {u!^ + Av'^k^ 


(5.9) 


d0i(xi)(502(a:2) 

Using the chain rule, it can be rewritten in the form 
d^Fniiw) 


(5.13) 


U 12 = 


dwidw2 


d{x — xi)S{x — X2), (5-14) 


Contraction of the vector indices in (15.81) leads to the 
following integrals over the momentum: 


/ 7:^Tw^4"ii(k)Pl8(k) (5.10) 

{27rYkd+v 12 V ; 781 ; V ^ 

I 

d^F^'^’^Ydwidw2 = 2s{wy-Y\wy [5i2w‘^ + 

+ 2ls{w^y~^ {\wy~^ {wi\2 


where = WkWk and (Xw) = XkWk- 

Now we have to contract the vertex factor (I5.16|) with 
the expression ()5.12l) . In order to find the renormaliza¬ 
tion constant, it is sufficient to retain only the terms pro¬ 
portional to the principal monomial in (15.5p and discard 
all the terms containing the factors of = A^A^. Com¬ 
bining all the relevant factors finally gives 


where 




(5.15) 


with the subsequent substitution Wi 9i{x). 
The differentiation gives 


2(s - l)uiiiU 2 ] + — l){w'^Y{XwY ^AiA 2 -I- 

+ W 2 A 1 ), 


(5.16) 


and g is defined in p.l5l) . 

The requirement that the renormalized analog of the 
function (15.31) be UV finite in the MS scheme gives: 


Znl = <l- 


9 


{Qi + aQ2)\, (5.19) 


r„i(x) = Fnl{x) I 1 - 


9_Y 9 {Qi + aQi) 
.TO/ Syd{dF2) 


,(5.17) 


8 yd{d + 2) 

and the anomalous dimension jni = P/i In Zni is 


where 


'yni — 


9 


8 d{d + 2) 


{Qi + aQ2) 


(5.20) 


Qi — — Ti^n -|- d){d — 1) -|- (d -f 1) ^(Z -I- d -I- 2), 

Q2 = —n{n + nd — d)(d — 1) -I- Z(Z -I- d -I- 2) (5.18) 


(we recall that we already set u = w = 1). 

Using the general expression (14.111) . for the critical di- 
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mension at the fixed point (14.81) we obtain 


“degree of anisotropy”): 


^nl — nAg + 7*; — + 


y{Qi + aQ2) 
6 (d- l)(d + 2) 


(5.21) 


with Ag from (14.121) . In particular, for the scalar operator 
one arrives at the expression 


which is negative and decreases as a grows: 

dAno/da < 0. (5.23) 

As we will see in the next section, this means that the 
anomalous scaling is indeed present in our model and be¬ 
comes more strongly pronounced as the degree of com¬ 
pressibility increases. 

For a fixed n, the dimensions (|5.2ip exhibit a kind of 
hierarchy with respect to the rank I (which measures the 

I 


dAni/dl > 0. (5.24) 

In contrast to the Gaussian model (see, e.g., [2i,il), 
this hierarchy becomes more strongly pronounced as a 
increases: 

d'^Ani/dWa > 0. (5.25) 

VI. OPERATOR PRODUCT EXPANSION AND 
THE ANOMALOUS SCALING 

A. General discussion and isotropic case 

The quantities of interest are, in particular, the pair 
correlation functions of the (renormalized analogs of 
the) operators (|5.1I) . In the following, we restrict our¬ 
selves with the equal-time correlations because they are 
Galilean invariant and do not bear strong dependence on 
the IR scale m = L~^ caused by the sweeping of small- 
scale vortices by the large-scale ones. Then one can write 


= n‘^‘"v‘^Pr]ni^qj{y,r,mr,c/yv) ~ (nr) 




( 6 . 1 ) 


where r = |r| = |x' — x|. The first equality follows from 
simple dimensionality considerations; then d‘f< and dp 
are the canonical dimensions of the correlation function, 
given by simple sums of the corresponding dimensions 
of the operators, and r]{- ■ ■) is some function of com¬ 
pletely dimensionless parameters. We have expressed the 
right-hand side in renormalized variables, when the ref¬ 
erence mass ^ is the substitute of the UV momentum 
scale A. The second (approximate) equality is valid in 
the IR asymptotic range /rr ^ 1 and follows from solv¬ 
ing the RG equation in the presence of the IR attrac¬ 
tive fixed point; A„j and Aqj being the critical dimen¬ 
sions of the operators from the left-hand side, given by 
Eq. (|5.21|) . In the following, we omit the RG invariant 
variable c(r) = c(/rr)^“/(/r^), which is restricted in the 
IR range; for more explanations, see [i^. We will also 
omit the indices of the scaling functions rj and C and do 
not display the dependence on the parameters v, y. that 
are fixed for a given physical setup. 

The inertial range corresponds to the additional in¬ 
equality mr <C 1. The form of the functions C in (EU) is 
not determined by the RG equations alone; they should 
be augmented by the OPE. In the case at hand, the 
OPE states that the equal-time product Fni{x)Fqj{x') 
at X = (x -I- x')/2 = const and r = x' — x —>■ 0 can be 
represented in the form 

Fr,i{x)Fqj{x') ^CF(r)F(i,x). (6.2) 

F 


Here the functions Cp are the Wilson coefficients, regular 
in m^, and F are, in general, all possible renormalized 
local composite operators allowed by the symmetry of the 
model and of the left-hand side. In the case at hand this 
implies that only Galilean invariant operators contribute. 
If these operators have additional vector indices, they are 
contracted with the corresponding (additional) indices of 
the coefficient functions Cp. 

Without loss of generality, it can always be assumed 
that the expansion (16.2|) is made in the irreducible op¬ 
erators with definite critical dimensions Ap. The cor¬ 
relation functions (lO) are obtained by averaging the 
expression ()6.2I) with the weight exp4S($), where 5(4>) 
is the (renormalized) action functional (13.5p . Then the 
quantities (F) appear on the right-hand sides. Consider 
first the isotropic case, then only the contributions from 
scalar operators survive. Their asymptotic behavior for 
m —>■ 0 is found from the RG equations for the operators 
F and has the form (F) oc (we recall that A^ = !)■ 

Thus, combining the expressions dSI]) and (16.2p gives 
the following inertial-range asymptotic representation for 
the scaling functions C: 

({mr) ~ E Ap{mr)^^, (6.3) 

F 

where all the coefficients Ap = Ap{mr) are regular in 
(mr)'^. 

Singularities for mr —^ 0 (and thus the anomalous scal¬ 
ing) result from the contributions in (16.31) of the operators 
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with negative critical dimensions, termed “dangerous” in 


41| . Clearly, if the number of such operators were fi¬ 


nite, the leading contribution would be determined by 
the operator with the lowest dimension. However, one 
can argue that, if at least one dangerous operator ex¬ 
ists in a model, the latter necessarily involves an infinite 
set of dangerous operators, and the spectrum of their di¬ 
mensions is not bounded from below; see Appendix A 
for discussion. In our case, from the expression (15.221) we 
can see that all the scalar operators Fno are dangerous, 
and their dimensions A„o increase without bound as n 
grows. 

Fortunately, the linearity of the original equation (12.61) 
in the field 9 imposes the restriction that the number of 
the fields 9 in all the composite operators in the expan¬ 
sion (lO) cannot exceed their number in the left-hand 
side; cf. the remark below Eq. (IQ) in Sec. El In turn, 
this means that only a finite number of the operators of 
the type can contribute to any given OPE. For the 
product (16.21) . these are the operators with k < n q. 
Thus, 


n+q 

^{mr) ~ ^ Ak{mr) + ... (6.4) 

k=0 

with Afco from (15.221) : the ellipsis stands for the “more 
distant” corrections to the small-mr behavior, given by 
the operators with derivatives and other types of fields. 
The leading term in (16.41) is determined by the operator 
with the maximum possible k = n + so that the final 
leading-order asymptotic expression for the correlation 
function dSH) in the inertial range fir ^ mr <C 1 has 
the form 


{FniFqj) {fir) (6.5) 

As already mentioned in the preceding section, the in¬ 
equality (15.231) means that the anomalous scaling be¬ 
comes more strongly pronounced as the degree of com¬ 
pressibility grows. We also note that the inequality 

An-t-iy ^ A^ -)- Ag, 

which follows from the explicit expressions (15.211) and, in 
fact, is required by the probabilistic theory, shows that 
the expression (16.5|) diverges for r —>■ 0. 


Uij ... Ui, -I- ..., where the ellipsis stands for the appro¬ 
priate subtractions with the Kronecker <5 symbols that 
make it irreducible. Upon substitution into the OPE 
for the product of two scalar operators, their tensor in¬ 
dices are contracted with the corresponding indices of 
the coefficient functions CF{r). This gives rise to the 
Gegenbauer polynomials, the d-dimensional analogs of 
the Legendre polynomials Pi (cos d), where d is the angle 
between the vectors r and n. For general anisotropy, the 
d-dimensional spherical harmonics appear on the right- 
hand side, while for products of tensor operators, their 
tensor analogs arise. 

Consider, as the simplest example, the pair correla¬ 
tion function (ISTl) of two scalar operators in the inertial 
range: 


(F„o(i,x)Fqo(t,x')) r x 

N 

X ^ c;Pi(cos'd)(mr)^^'-I-... (6.6) 

with N = n + q and Ajvz from (I5.2ip ; ci are numerical 
coefficients and the ellipsis stands for the “distant” con¬ 
tributions with I > N. The inequality (15.241) means that 
the anisotropic contributions in (16.61) exhibit a kind of 
hierarchy related to the degree of anisotropy 1 : the lead¬ 
ing contribution is given by the isotropic “shell” {I = 0), 
while the contributions with I > 1 give only corrections 
which become relatively weaker as mr —> 0, the faster the 
higher the degree of anisotropy I is. Similar hierarchy, 
observed earlier in numerous models, e.g., [l9l - l^ . [SJ- 
IHiIHIj gives quantitative support for Kolmogorov’s hy¬ 
pothesis of the local isotropy restoration. 

The inequality (15.251) means that the hierarchy (15.241) 
becomes stronger as the degree of compressibility a 
grows: the anisotropic corrections are getting further 
from one another and from the isotropic term, contrary to 
the situation observed earlier for passive vector field, ad- 
vected by Kraichnan’s ensemble [22j| . A similar discrep- 
ancy for the scalar field was encountered recently in |48l | . 
This means that the results obtained on the base of sim¬ 
plified “synthetic” ensembles should be taken with some 
precaution. 


VII. CONCLUSION 


B. Effects of large-scale anisotropy 

Consider effects of the anisotropy, introduced into the 
system at large scales ^ L through, say, the large-scale 
field = UiB^ or through the correlation function of the 
artificial random force. Then the irreducible tensor com¬ 
posite operators acquire nonzero mean values, built of 
the vector n: for example, the mean value of the second- 
rank operator is proportional to the irreducible tensor 
UiUj — 5ijld. In general, the mean value of any 1-th 
rank irreducible operator is proportional to the tensor 


We have studied the model of passive vector field, 
advected by a turbulent flow. The latter is described 
by the Navier-Stokes equations for a strongly compress¬ 
ible fluid (12.11) . (12.21) with an external stirring force with 
the correlation function oc k'^' (HH), (113]). From 

physics viewpoints, the model describes magnetohydro¬ 
dynamic turbulence in the so-called kinematic approxi¬ 
mation, where the effects of the magnetic field on the 
dynamics of the fluid are neglected. 

The full stochastic problem can be cast as a field the¬ 
oretic model with the action functional specified in (12.81) 



















12 


and That model appears multiplicatively renor- 

malizable, so that the corresponding RG equations can 
be derived in a usual way. They have the only IR attrac¬ 
tive fixed point in the physical range of parameters, so 
that various correlation functions reveal scaling behavior 
in the IR region. 

Their inertial-range behavior was studied by means of 
the OPE; existence of anomalous scaling (singular power¬ 
like dependence on the integral scale L) was established. 
The corresponding anomalous exponents were identified 
with the scaling (critical) dimensions of certain compos¬ 
ite fields (composite operators), namely, powers of the 
magnetic field. They can be systematically calculated as 
series in the exponent y. The practical calculation was 
accomplished in the leading order; the results are pre¬ 
sented in (j5.21|) . 

The results obtained are quite similar to those derived 
earlier for the vector fields advected by synthetic veloc¬ 
ity ensembles [l^, [2^: the anomalous scaling becomes 
more remarkable as the degree of compressibility a in¬ 
creases; the anisotropic contributions form an hierarchy 
related to the degree of anisotropy I, so that the leading 
inertial-range contribution is the same as for the isotropic 
case. However, that hierarchy becomes stronger as the 
degree of compressibility grows, in contrast to what was 
observed in [22j| for the Kraichnan’s rapid-change ensem¬ 
ble. In this respect, our results are close to what was 
recently observed for the scalar field, advected by the 
same velocity ensemble (d^ . 

From a more theoretical point of view, it is important 
that in our case, the anomalous exponents are associated 
with the critical dimensions of certain individual com¬ 
posite operators, exactly as in the RG-I-OPE treatment 
of the rapid-change models; see, e.g., [38l - l40l | . In 

the zero-mode approach to the latter, the anomalous ex¬ 
ponents are related to the so-called zero modes (unforced 
solutions) of the exact differential equations satisfie d by 
the equal-time correlation functions; see, e.g., [iTl - llftl^ . 
In a more general sense, zero modes can be interpreted 
as certain statistical conservation laws in the dynamics 
of particle clusters [s^. The close resemblance in the 
RG-bOPE pictures of the origin of anomalous scaling for 
the present model and its rapid-change predecessors sug¬ 
gests that the concept of zero modes (and thus that of 
statistical conservation laws) is also applicable in much 
more realistic models. 


Appendix A: Infinite number of negative dimensions 

Let E'(x) be a certain renormalized composite opera¬ 
tor in a certain multiplicatively renormalizable field the¬ 
oretic model with an IR attractive fixed point of the RG 
equations. Assume that F has a definite negative critical 
dimension, < 0, and assume that it is the lowest di¬ 
mension in the model (that is, F is the “most dangerous” 
operator). 


Gonsider its pair correlation function: 

/ dk 

^^^£)i.(fc)exp{ik-(x'-x)}, (Al) 

where k = |k|. Our assumption implies that in the IR 
range the function Dp(k) has the asymptotic form 

DF[k) ~ /c-'^+2+^^/(m/A:), (A2) 

where m is some typical IR momentum scale, f{m/k) 
is some dimensionless scaling function, and we assumed 
that Am = 1 (like in our model). Now consider the mean 
value 

r dk 

The function f(rn/k) provides IR regularization of the 
integral (IA3I) . The question is whether this integral re¬ 
mains convergent for large k if the exact (unknown) func¬ 
tion Dpik) is replaced with its asymptotic form (IA2I) . 
According to the OPE, the asymptotic behavior of the 
function f(rn/k) for large k, or, equivalently, for small m 
is determined by the contribution of the most dangerous 
operator, which by assumption is F itself: 


f{m/k) ~ {m/k)^’". 

(A4) 

Thus, for large k we have 


Dpik) ~ 

(AS) 


and the integral (IA3I) remains convergent upon the sub¬ 
stitution of (IA2I) . In turn, this means that it is expressed 
only in terms of the IR parameter m (UV regularization 
A can be removed). Then it is easily found from the 
dimension: 

(F2(x)) ~ (A6 ) 

Expression (IA6|) means, however, that the operator 
has the negative dimension 2 IS.f < Ap < 0 which is 
smaller than that of F. We arrive at the contradiction 
with our initial assumption about F. 

To avoid possible misunderstanding we stress that our 
consideration does not mean that F^ is necessarily dan¬ 
gerous and its dimension is exactly 2Ai7’ (although this 
indeed happens, e.g., for the powers of the velocity field 
in the stochastic NS problem; see [4l|). It means that 
operators with negative dimensions, if any, always ap¬ 
pear in a model as infinite families, with the spectrum 
of dimensions not bounded from below. This fact should 
be taken into account in axiomatic or phenomenolo gica l 
implementations of the OPE to models of turbulence . 

ACKNOWLEDGMENTS 

The authors are indebted to L. Ts. Adzhemyan, 
Michal Hnatich, Juha Honkonen, and M. Yu. Nalimov 
for discussion. 






13 


The Authors acknowledge the Saint Petersburg State Mariia Kostenko was also supported by the Dmitry 
University for the research grant 11.38.185.2014. Zimin’s Dynasty Foundation. 


[1] U. Frisch, Turbulence: The Legacy of A. N. Kolmogorov 
(Cambridge University Press, Cambridge, UK, 1995). 

[2] P. S. Iroshnikov, Sov. Astron. 7, 566 (1964). 

[3] R. H. Kraichnan, Phys. Fluids 8, 1385 (1965). 

[4] L. F. Burlaga, J. Geophys. Res. 96, 5847 (1991); 97, 4283 
(1992); 

E. Marsch and S. Liu, Ann. Geophys. 11, 227 (1993). 

[5] R. Grauer, J. Krug, and C. Marliani, Phys. Lett. A 195, 
335 (1994). 

[6] G. Einaudi, M. Velli, H. Politano, and A. Pouquet, As- 
trophys. Journ. 457, L113 (1996). 

[7] R. Grauer and C. Marliani, Phys. Plasmas 2, 41 (1995); 
Physica Scripta T 67, 38 (1996). 

[8] G.-Y. Tu and E. Marsch, Space Sci. Res. 73, 1 (1995). 

[9] G. Pagel and A. Balogh, Nonlin. Processes in Geophysics 
8, 313 (2001). 

[10] R. Bruno, V. Carbone, B. Bavassano, L. Sorriso-Valvo, 
and E. Pietropaolo, Mem. Sos. Astrophys. It. 74, 725 
(2003). 

[11] R. Bruno, B. Bavassano, R. D’Amicis, V. Carbone, 
L. Sorriso-Valvo, and A. Noullez, Geophys. Research Ab¬ 
stracts 9, 08623 (2003). 

[12] G. Salem, A. Mangeney, S. D. Bale, and P. Veltri, Astro¬ 
phys. J. 702, 537 (2009). 

[13] P. D. Mininni and A. Pouquet, Phys. Rev. E 80, 025401 
(2009). 

[14] M. W. Kunz et ah, arXiv:1501.06771 astro-ph.HE]. Sub¬ 
mitted to J. Plasma. Phys. (2015). 

[15] L. Sorriso-Valvo et ah, arXiv:1505.97879 physics.space- 
ph]. 

[16] C. C. Lalescu et ah, Phys. Rev. Lett. 115, 025001 (2015). 

[17] M. Vergassola, Phys. Rev. E 53, R3021 (1996). 

[18] 1. Rogachevskii and N. Kleeorin, Phys. Rev E 56, 417 
(1997). 

[19] A. Lanotte and A. Mazzino, Phys. Rev. E 60, R3483 
(1999). 

[20] 1. Arad, L. Biferale, and 1. Procaccia, Phys. Rev. E 61, 
2654 (2000). 

[21] N. V. Antonov, A. Lanotte, and A. Mazzino, Phys. Rev. 
E 61, 6586 (2000). 

[22] N. V. Antonov, J. Honkonen, A. Mazzino, and P. Mura- 
tore Ginanneschi, Phys. Rev. E 62, R5891 (2000). 

[23] L. Ts. Adzhemyan, N. V. Antonov, A. Mazzino, P. Mu- 
ratore Ginanneschi, and A. V. Runov, Europhys. Lett. 
55, 801 (2001); 

H. Arponen, Phys. Rev. E, 79, 056303 (2009). 

[24] N. V. Antonov, M. Hnatich, J. Honkonen, and M. 
Jurcisin, Phys. Rev. E 68, 046306 (2003). 

[25] N. V. Antonov and N. M. Gulitskiy, Phys. Rev. E 85, 
065301(R) (2012); E 87, 039902(E) (2013); 

Phys. Rev. E 91, 013002 (2015); 
arXiv: 1506.05615 (2015). 

[26] E. Jurcisinova and M. Jurcisin, J. Phys. A: Math. Theor. 
45, 485501 (2012). 

[27] E. Jurcisinova and M. Jurcisin, Phys. Rev. E 88, 011004 
(2013) 

[28] E. Jurcisinova, M. Jurcisin, and Remecky, Phys. Rev. E 


88, 011002 (2013); 

E. Jurcisinova, M. Jurcisin, and P. Zalom, Phys. Rev. E 

89, 043023 (2014). 

[29] N. V. Antonov and N. M. Gulitskiy, Theor. Math. Phys. 
176(1), 851 (2013). 

[30] J. D. Fournier, P. L. Sulem, and A. Pouquet, J. Phys. A 
15, 1393 (1982); 

L. Ts. Adzhemyan, A. N. Vasil’ev and M. Gnatich, Theor. 
Math. Phys. 64, 777 (1985); 72, 940 (1987). 

[31] A. P. Kazantsev, Sov. Phys. JETP 26, 1031 (1968). 

[32] R. H. Kraichnan, Phys. Fluids 11, 945 (1968); Phys. Rev. 
Lett. 72, 1016 (1994). 

[33] G. Falkovich, K. Gaw§dzki, and M. Vergassola, Rev. 
Mod. Phys. 73, 913 (2001). 

[34] N. V. Antonov, Phys. Rev. E 60, 6691 (1999); 

L. Ts. Adzhemyan, N. V. Antonov, and J. Honkonen, 
Phys. Rev. E 66, 036313 (2002). 

[35] N. V. Antonov, Physica D 144, 370 (2000). 

[36] M. Hnatich, E. Jurcisinova, M. Jurcisin, and M. Repasan, 
J. Phys. A: Math. Gen. 39, 8007 (2006). 

[37] M. Holzer and E. D. Siggia, Phys. Fluids 6, 1820 (1994). 

[38] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil’ev, 
Phys. Rev. E 58, 1823 (1998) 

Theor. Math. Phys. 120, 1074 (1999). 

[39] L. Ts. Adzhemyan, N. V. Antonov, V. A. Barinov, Yu. S. 
Kabrits, and A. N. Vasil’ev, Phys. Rev. E 63, 025303(R) 
(2001); E 64, 019901(E) (2001); E 64, 056306 (2001). 

[40] N. V. Antonov, J. Phys. A: Math. Gen. 39, 7825 (2006). 

[41] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil’ev, 
Sov. Phys. JETP 68, 733 (1989); 

L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil’ev, 
Usp. Fiz. Nauk, 166, 1257 (1996) [In Russian, Engl. 
TransL: Phys.-Usp., 39, 1193 (1996)]; 

L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasiliev, 
The Field Theoretic Renormalization Group in Fully De¬ 
veloped Turbulence (Gordon & Breach, London, UK, 
1999). 

[42] S. S. Moiseev, A. V. Tur, and V. V. Yanovskii, Sov. Phys. 
JETP 44, 556 (1976). 

[43] V. S. L’vov and A. V. Mikhailov, Sov. Phys. JETP 47, 
756 (1978). 

[44] 1. Staroselsky, V. Yakhot, S. Kida, and S. A. Orszag, 
Phys. Rev. Lett. 65, 171 (1990). 

[45] L. Ts. Adzhemyan, M. Yu. Nalimov, and 

M. M. Stepanova, Theor. Math. Phys. 104, 305 
(1997). 

[46] D. Yu. Volchenkov and M. Yu. Nalimov, Theor. Math. 
Phys. 106, 375 (1996). 

[47] N. V. Antonov, M. Yu. Nalimov, and A. A. Udalov, 
Theor. Math. Phys. 110, 305 (1997). 

[48] N. V. Antonov and M. M. Kostenko, Phys. Rev. E 90, 
063016 (2014). 

[49] L. D. Landau and E. M. Lifshitz, Fluid Mechanics 2nd 
ed. (Pergamon Press, Oxford, UK, 1987). 

[50] H. K. Moffat, Magnetic field generation in electrically 
conducting fluids (Cambridge University Press, Cam¬ 
bridge, 1978); 



14 


Yu. V. Novozhilov and Yu. A. Yappa, Electrodynamics namics (Boca Raton, Chapman & Hall/CRC, FL, 2004). 

(Mir Publishers, Moscow, 1981). [53] A. M. Polyakov, Nucl. Phys. B 396, 367 (1993); 

[51] J. Zinn-Justin, Quantum Field Theory and Critical Phe- G. Falkovich and A. Zamolodcikov, J. Phys. A: Math. 

nomena (Clarendon, Oxford, UK, 1989). Theor. 48, 18FT02 (2015). 

[52] A. N. Vasiliev, The Field Theoretic Renormalization 
Group in Critical Behavior Theory and Stochastic Dy- 



